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Abstract 

A numerical study of the Faddeev equation for bosons is made with two-body interactions at or close 
to the Unitary limit. Separable interactions are obtained from phase-shifts defined by scattering length and 
effective range. In EFT-language this would correspond to NLO. Both ground and Efimov state energies are 
calculated. For effective ranges ro > and rank-1 potentials the total energy Et is found to converge with 
momentum cut-off A for A >'~ 10/ro . In the Unitary limit (1/a = ro = 0) the energy does however diverge. 
It is shown (analytically) that in this case Et ~ Eut^ ■ Calculations give Eu = —0.108 for the ground state 
and i?u = —1. X lO"** for the single Efimov state found. The cut-off divergence is remedied by modifying 
the off-shell t-matrix by replacing the rank-1 by a rank-2 phase-shift equivalent potential. This is somewhat 
similar to the counterterm method suggested by Bedaque et al. This investigation is exploratory and does 
not refer to any specific physical system. 
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Systems involving particles with 2-body interactions at or close to the Unitary limit have become of specific 
interest in the physics community during the last several years for reasons that have repeatedly been pointed out 
in numerous publications related to both atomic and sub-atomic problems. The Unitary limit is here defined 
as being that for which the scattering length a and effective range ro are infinite and zero respectively. In that 
case the only scale left in a system of fermions of 'infinite' extension vifould be the Fermi-momentum. The total 
energy would then have to be proportional to the Fermi-energy as that is the only scale left in this problem. 
Experimental results point to Et OAAEp with Ep being the kinetic energy of the 'unperturbed' kinetic 
^ , energy of the system. Theoretical determination of the this universal constant is a many-body problem. It 

(— I ' should however only involve some simple constants and might apriori seem straightforward. It has however 

been found to be a theoretical challenge. Some calculations of the author using Brueckner methods show for 
example a very strong dependence on the effective range rg. There is also in any theoretial calculation a 
necessary cut-off A in momentum-space that renders the 2-body interaction a function of this A. This does 
\^ ' introduce another length parameter into the theory. In the previous work of the author, and shown below, the 

^SJ , rank-1 separable interaction is in the Unitary limit singular at the momentum A. This fact introduces another 

0^ ' problem; the applicability of a many-body theory in this limit. The applicability of the Bruckner method was 

for example questioned after the realisation of the very large correlations and consequently, "wound-integral" 
04 ' in this case. The "model" space represented by a zero-temperature Fermi distribution is no longer adequate. A 

Green's function method including spectral broadening would be more appropriate. 

As opposed to the 'infinite' system the three-body system is exactly solvable by the Faddeev method [1]. It 
therefore presents a more interesting project for a numerical study with interactions at and near the Unitary 
limit. Of specific interest here is also the phenomena first brought to the attention by Vitali Efimov. % '7 He 
^> , found the surprising fact that bosons interacting with a resonance in the 2-body state (i.e. 1/a 0) would 

^ ' result in a strongly bound three-body system and with a spectrum of loosely bound excited states. 

5^ I The inverse scattering method as applied here uses two-body on-shell data (scattering length and effective 

range) as input. Although on-shell properties of the t-matrix are fitted exactly, many-body calculations involve 
also off-shell t-matrix elements. These are not derivable from experimental two-body data. Efeo even if a rank- 
1 potential is sufficient to fit the on-shell data as is often the case it leaves the off-shell undetermined. The 
extension to a rank-2 will allow for a phase-shift equivalent interaction with different off-shell t-matrix elements. 
This provides a practical tool for exploring off-shell or equivalently, three-body effects. This method will be 
used in some cases below. It was used by the author in earlier work, (see e.g. |TT]) The on-shell data relate to 
the asymptotic form of the two-body scattering wave-function. These have to be preserved when increasing the 
rank which implies that the interaction should be modified at short range in coordinate space as shown below. 

In the present report we focus only on the ground and Efimov state energies, as well as on the dependence on 
scattering lengths and effective ranges and on the questions of convergence as a function of cut-off in momentum- 
space. 
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^This presents a problem, of course not restricted to the use of the inverse scattering method. 
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In Section 2 is found a presentation of the necessary tools which are the Faddeev equation and the inverse 
scattering method. Section 3 show resuhs of numerical calculations in 3 subsections wth 8 Figures. Section 4 
is a summary and some discussion of the results. 



2 Formalism 

The Faddeev three-boson equation for a spin-independent rank-1 separable attractive potential V{k, k') = 
—v{k)v{k') is given by 

2 /-^ ^(|k+lq|M|q-|-lk|) ^ 



with 
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Z(s) = l + -^/ v^{k){s-k^)-^k^dk (2) 



The extension to the rank-2 potentials used below in Section 3.2 is straight-forward. [4", "5^ 
The separable interactions are calculated from phase-shifts by an inverse scattering method that dates back 
at least some 40 years. [5] Some recent applications by the author can be found in refs where details can be 
found. [71 m [9l [Ml [ini E] where some details of the method are also shown. The input are phase-shifts which 
in general can be either experimental or otherwise defined. For the purpose of this investigation they will be 
defined by a scattering length a and and an effective range vq. A rank-1 separable potential is then sufficient 
to reproduce the input phase-shifts exactly. (If the phase-shift were to changes sign such as in the nuclear ^Sq 
case a rank two potential is necessary. jTj). As mentioned above a higher rank may be required in order to 
accomodate three-body data. The present work is not specific to any particular system other than that the 
2-body interaction is close to the Unitary limit for which universality would apply. It is shown however that an 
off-shell correction (three-body force) via a rank-2 potential can be used to prevent the ultraviolet divergence 
and collapse of the three-boson system in that limit. 
In the case of a rank-1 attractive potential one has 

V{k,k') ^ ^v{k)v{k') (3) 

Inverse scattering then yields (e.g. ref [7l[6]) 

v^{k) = ^^sinS{k)\D{k^)\ (4) 



where 



k 



2 ^ 

tt' ./o k^-k'^ 



(5) 



where V denotes the principal value and 6{k) are the phaseshifts. A provides a cut-off in momentum-space and 
the interaction is fully defined by the phase-shifts, the two-body binding energy Eb and by A. The effect of the 
cut-off will be exploited below. Eb is calculated from 

^^iEB) = - + ItoEb (6) 
a 2 

It is set to zero for unbound states. 

For the rank-2 potentials that are used in Sect. 3.2 the method of Chadan and Fiedeldey was applied. [T^ [12] 
(see also ref. |7]). In this a set of initial phaseshifts is assumed to be provided. An arbitrary interaction Vi, 
is then assumed, defined either by another set of phases or explicitly. A second potential, V2 can then be 
constructed so that the rank-2 potential given by Vi and V2 reproduce the initial phases. If the initial phases 
are, as in our case below, all of the same sign, they can be reproduced by a rank-1 potential from eq. ([4]). The 
extension to a rank-2 potential, allows for an arbitrary change of off-shell behaviour or in other words of the 
three-body term while preserving the fit to the initial two-body phase-shifts. 

With S{k) — ^, the unitary limit, these phases can be reproduced by a rank one-potential and one finds 
(e.g.lH) 
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where the factor A (= 1 in the unitary hmit) is introduced for later presentation of results where the three-boson 
binding is calculated as a function of this strength factor. Note that ^^(fc) — oo for fc — A. 
Note also that for A fc one finds 

vlik) ^ (8) 

In this limit, but only in this limit , the unitary interaction will then be a (5- function in coordinate space with 
the strength inversely proportional to the cut-off. But eq. ([7]) shows that a finite A results in an abrupt increase 
in strength and a singularity at A: = A to preserve the condition 5 = ^ for all fc < A. 

With the interaction ^ and with the momenta in units of the cut-off A {k\ — k/A and qa = q/A) the 
Faddeev equation is: 

2A /"^ w(|kA + 5qA|)t^(|qA + 5kA| 



^^'^ ^ AE^W - hi) Jo .i + qA-kA + Ai-K(A) ^^'^^^^^ 

where Eu{X) = Et/A^- 

The function I{s) is after a change of variables [/ca — sin{0)\ given by 

2A 

I{s) = 1 + —V / sin^{e){s - sin^{e))-^de (10) 

Jo 

Note that any interaction for which w^(fc) = F(k)/A would result in a similar universal equation with 
Et oc A2. 

It is also important to note that with s = fc^ and A = 1 one finds 

I{s) = 

for ALL < s < 1. This leads to a Reactance matrix element 

< klTZlk >= ytanSik) oo 
k 

i.e. 5{k) = ^ for < fc < A which is the condition for a Unitary interaction with a cut-off A in momentum-space. 

The associated resonance is of course also the origin og the Efimov-states in the three-body system. The 
number of such states were predicted to be [2] 

iV = (l/7r)Zn(|a|/ro) (11) 

The problem with the renormalisation of the three-boson system (i.e. the A^ divergence) was addressed 
by Bedaque, Hammer and van Kolck.^16; They resolve it by introducing a three-body counter-term. Here 
this is done by extending the rank-1 potential to a rank-2 by the method described above. The 'arbitrary' 
interaction Vi is defined by repulsive phases Sk = —Tc * k, k being the relative momentum and Vc a constant 
determined below. This simulates a short-ranged repulsion which effectively removes the ultra-violet divergence 
experienced with the rank-1 potential. The off-shell t-matrix elements in the Faddeev equation are affected by 
this modification of the interaction with results shown below. 



3 Numerical Results 

The energy of the three-boson system was calculated by solving the Faddeev equation numerically either by 
iteration or by matrix diagonalisation in the conventional fashion. The separable interaction was obtained by 
the inverse scattering method referred to above for a range of scattering lengths and effective ranges including 
the Unitary limit, In the results presented below the energies are in units of jm and the lengths in units of 
/m, but are in general universal. 

3.1 Ground States with rank-1 potential 

It was verified numerically that the quantity £'„(A), defined above, is independent of the cut-off A. The 
dependence of A for the ground state is shown in Fig. [T]with ^^^(A = 1) = —0.108 (i.e. the Unitary limit) while 
E^ = for A ~ 0.77 with a nearly linear dependence on A. 
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Figure 1: The energy Eu, defined in the text, as a function of the strength A of the Unitary interaction ([Tj- The 
energy of the three-boson system would be Et = —Eu{X) x A^. The calculated values are indicated by points 
and connected by lines for clarity. 




Figure 2: The energy of the three-boson system is shown as a function of the scattering length and three effective 
ranges that are the parameters defining the 2-body interaction as described in the text. The calculated values 
are indicated by points and connected by lines for clarity. 
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Figure 3: The energy Et of the three-boson system is shown as a function of momentum cut-off A for 
three different values of scattering lenghts and four different effective ranges as indicated. See text for further 
discussion. The calculated values are indicated by points and connected by lines for clarity. 



Fig. [5] shows the energy as a function of 1/a for three chosen values of effective ranges tq. One notices a 
very strong dependence on the effective range. 

Fig. |3] shows energies Et as a function of cut-off A for 1/a = and some selected values of effective ranges 
tq. It was shown above that for 1/a — and tq = 0, i.e. the Unitary limit , Et — O.IOSA^. This is shown 
numerically by the curve labelled " .0" . The figure shows that the convergence with A improves as the effective- 
range is increased. The value of A = Ac at which the asymptotic value is reached scales roughly as 10/ro. 

The size of the three-boson system in momentum-space scales roughly the same, consistent with that the 
size in coordinate space would be w tq. Fig. |4] shows the rms radius Rrms of the ^(q) function. It follows 

closely the A-dependence of Et shown in Fig. [3l H One concludes that the size of the system (in momentum- 
space) determines the minimum range Ac of momenta that the interaction has to span. This is analogous to the 
situation found in nuclear matter Brueckner calculations where the minimum cut-off is found to be 2k f i.e. twice 
the fermi-momemtum. [9; In the present work on the three-boson system we also find, quite naturally, that Rrms 
is inversely proportinal to ro. This is a difference from the Brueckner calculations where the effect of correlations 
on the momentum distribution is ignored and approximated by the non-interacting fermi-distribution. This is 
an approximation related to the quasi-particle approximation which is implicit in the Brueckner method. The 
Green's function approach goes beyond this approximation including the finite width of the spectral-functions 
and the momentum-distribution is then wider. 

3.2 Rank-2 potential in Unitary Limit 

The renormalisation of the non-relativstic three-body problem with short ranged forces was addressed by Be- 
daque et al[T2]. As shown above, the three-boson system with the two-body unitary interaction (O collapses 
as A^. Bedaque et al suggests to introduce a three-body counterterm for the renormalisation. One may al- 
ternatively choose to introduce a similar counterterm by changing the off-shell two-body t-matrix. As alreday 
announced in Sect 2 this is done by replacing the rank-1 by a rank-2 interaction with Vi chosen so as to modify 
the short-ranged, ultraviolet, part of the interaction in order to prevent the collapse. With 5k = — Tc * k the 

^One may however note a slower approach to the asymptotic value at large A. This is a general characteristic of any energy vs 
size display. 
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Figure 4: The 'size' of the three-boson system in momentum-space is shown as a function of cut-off A for the 
same effective ranges as in Fig. [3] Only the results for 1/a = are shown. Note the similar A-dependence as is 
exhibited in Fig. |3l The calculated values are indicated by points and connected by lines for clarity. 

related potential Vi is obtained by inverse scattering. Fig. [S] shows the three-boson energy as a function of 
cut-off A for three different values of Tc- The lowest curve is for Tc = i.e. the rank-1 potential and shows the 
A^ divergence. The upper five curves are for increasing values of Tc as shown in the figure caption. A drastic 
change is seen with Tc = 0.12 showing convergence for large A. 

3.3 Efimov States 

There are many publications related to the low-bound excited states found as solutions of the Faddeev equations, 
states first discoverd by and named after Vitali Efimov [5] Numerous calculations were done at and close to the 
unitary limit. Never was found more than one state that could be identified as an 'Efimov' state although eq. 
pT|) suggests several states should be found with 1/a — rg '~ 0. The search for these states was in general 
very elusive and in particular very sensitive to very small changes in the low momenta part of the interactions. 
The broken curve in Fig. [6] shows the excited state energy £'„(A) obtained from solving eq. ([9]). The solid line 
represents thw two-body bound state energy. The sole Efimov state is seen to coincide with the two-body at 
A = 1.08. The Faddeev equation also yields numerous three-boson energies located above the two-body curve. 
These represent dimers, two bound bosons, and a free boson. The ground state energies are some fifty times 
deeper. 

Another example is shown in Fig. [71 The separable interaction is here defined by a rank-1 potential with 
scattering length a = —2 and an effective range tq = .03. Solutions of the Faddeev equation are shown as 
a function of the strength multiplier A. The solid line shows the two-boson bound state energies. (Not clearly 
shown is that they are unbound for A < 1.04). The sole line below this is the Efimov state. The numerous lines 
above are dimer-|-l states. 

Yet another example of the numerous calculations that were made is shown in Fig. [8l The effective range 
is here rp =0.1 and the scattering length is allowed to vary over the range indicated. The dots connected with 
broken lines indicate solutions of Faddeev equations with the binding energy Eb from eq. (I6|). According to 
other works these should with incresing 1 /a approach the bound dimer line. This is not the case here. This may 
be a characteristic of the separable interaction. Humberston et al [17] show a comparison of separable (non- 
local) and local (Yukawa and exponential) interactions used in three-boson calculations. The energy increases 
much faster with the strength of the interaction for the separable than it does for the local interaction. In 
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Figure 5: The three-boson energy as a fuction of A for different values of r^. From bottom and up: = 
0, .001, .0012, .0013, .0014., 0015 [fm]. The calculated values are indicated by dots, connected by straight lines. 
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Figure 6: The broken line shows the Efimov state as a function of the strength A of the unitary interaction ([7]). 
The energy Eu is given in units of as in eq. ([9]). The full line shows the two-body bound state energy. 
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Figure 7: Energies as a function of the strength A. The sohd hue shows the two-boson bound state energies. 
The sole hne below this is the Efimov state. The numerous lines above are dimer+l states and would not be 
classified as Efimov states. 

order to investigate the effect of the binding, another set of calculations were made shown by the broken line 
connecting squares. Here Eb is set to zero in eq. ([5]). It does have a 1/a dependence similar to what was 
expected from the literature. Only one Efimov state is found here with rg — 0.1. From eq. pT|) 2 states could 
be expected only for 1/a < .025. It should be mentioned that expression PT|) was tested by Huber [TB]. Our 
result may be associated with the rank-1 potential. 

4 Summary and Comments 

Fig. [2] shows that the energy of the bound system of three-bosons depend strongly on the range tq of the 
two-body interaction. This is consistent with earlier results for the infinite fermion-system. [8l I15j Fig. |3] shows 
that as a function of the momentum cut-off A, the energies converge toward asymptotic values Et that, as in 
Figl2l are functions of the range .03 < tq < .1 but largely independent of the scattering length for \a\ ^ tq. 
Convergence was reached at Ac ~ 10/ro ^ 5 x Rrms- The size of the trinicr at equilibrium (saturation), the 
inverse of Rrms scales as ro and the range A has to scale with Rrms consistent with our results. 

For comparison the curve labelled by " .0" in Fig. [3] shows the energy vs A in the unitary limit. As shown 
above (eq. ©) this limit has to be treated as a special case giving the analytic result Et = Eu{X)A'^. The 
numerical result shown by Fig. [T]is £'„(A = 1) = -0.108 while Eu{X) = for A 0.77. 

The quadratic divergence and simaltaneous collapse of the boson trimer in the Unitary limit was dampend by 
a renormalisation procedure consisting of a counter-term represented by a rank-2 potential with results shown 
in Fig. [5l 

Efimov states were identified although not quite as expected which may be a consquence of the specific 
interactions. 
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